Electronic band structure and exchange coupling constants in ACr2X4 spinels 
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We present the results of band structure calculations for ACT2X4, {A—Zn, Cd, Hg and X=0, 
S, Se) spinels. Effective exchange coupling constants between Cr spins are determined by fitting 
the energy of spin spirals to a classical Heisenberg model. The calculations reproduce the change 
of the sign of the dominant nearest-neighbor exchange interaction Ji from antiferromagnetic in 
oxides to ferromagnetic in sulfides and selenides. It is verified that the ferromagnetic contribution 
to Ji is due to indirect hopping between Cr t2g and Cg states via X p states. Antiferromagnetic 
coupling between 3-rd Cr neighbors is found to be important in all the ACT2X4 spinels studied, 
whereas other interactions are much weaker. The results are compared to predictions based on the 
Goodenough-Kanamori rules of superexchange. 

PACS numbers: 71.20.-b, 71.70.Gm, 75.30.Et 



I. INTRODUCTION 

Chromium spinels provide unrivalled possibilities for 
studying magnetic interactions in solids. In these com- 
pounds with a general formula ACT2X4, where A is a 
divalent nonmagnetic cation (Mg, Zn, Cd, or Hg) and 
X is a divalent anion (O, S, or Se), a Cr'^+ ion is in 
the 3d^ configuration. Its three 3d electrons occupy the 
majority-spin states of a completely spin polarised t2g 
sub-shell leading to the total spin S=3/2. Although 
charge and orbital degrees of freedom in the ACT2X4 
spinels are frozen, these compounds show wide variety 
of magnetic properties ranging from those of a strongly 
frustrated antiferromagnet to a Heisenberg ferromagnet. 
Depending on the chemical composition their effective 
Curie- Weiss temperature (8cw) varies from —400 K in 
oxides to 200 K in selenides fii^ which indicates that the 
sign of the dominant exchange interaction changes from 
antiferromagnetic (AFM) to a ferromagnetic (FM) one. 

In the ylCr204 spinels AFM nearest-neighbor interac- 
tions between Cr spins residing on a pyrochlore lattice are 
geometrically frustrated. The magnetic ground state of 
a frustrated antiferromagnet is highly degenerate which 
leads to unusual low-temperature properties.'^'* Cr oxide 
spinels remain paramagnetic with the Curie- Weiss form 
of the magnetic susceptibility down to temperatures well 
below IGcwl of 398, 71, and 32 K for A=Zn, Cd, and 
Hg. ZnCr204 (TAr=12.5 K) and CdCr204 (Ta,=7.8 K) 
undergo a first order phase transition of the Spin-Peierls 
type into a magnetically ordered Neel state at tempera- 
tures much lower than the characteristic strength |Ocw| 
of the interaction between Cr spins. ^'^ The transitions 
are accompanied by cubic to tetragonal structural dis- 
tortions, however, the sign of the distortions and the 
magnetic order below Tn are different. In ZnCr204 the 
lattice contracts along the c axis (c < a) and the Neel 
state has a complex commensurate spin structure with 4 
characteristic wave vectors. In contrast, the lattice of 
CdCr204 expands below Tn (c > a) and its ordered 
state is incommensurate with a wave vector Q=(0,<5,1) 
with S ^0.09.^ HgCr204 also undergoes a transition 



to a magnetically ordered state at Tn=5.8 K but the 
symmetry of the lattice lowers to orthorhombic.^ Re- 
cently, a metamagnetic transition and a wide magnetisa- 
tion plateau with the magnetic moment equal to one-half 
of the full Cr moment have been observed in CdCr204 
and HgCr204P^'^'2^ 

In ACr2S(e)4 spinels dominant ferromagnetic inter- 
actions are not geometrically frustrated. Nevertheless, 
Cr spins in ZnCr2S4 and ZnCr2Se4 form helical spin 
structures below 15.5 and 18 K, respectively . ^ 1 1 ^ ^ 1 13 jjj 
ZnCr2S4 the helical structure coexists at low tempera- 
tures with a coUinear AFM one. The transitions into 
helically ordered state are supposed to occur because of 
competing FM nearest-neighbor interaction and AFM 
interactions between more distant Cr neighbors. Re- 
cently, it has been shown that CdCr2S4 and HgCr2S4 
exhibit ferroelectric behavior with strong increase of the 
dielectric constant below the temperature of magnetic 
ordering li^ii^ 

The diversity of magnetic properties of the ACT2X4 
spinels can hardly be explained without understand- 
ing the mechanism of exchange interactions between Cr 
spins, their range, and relative strengths. So far theoret- 
ical analyses of the effective exchange interactions in Cr 
spinels were mostly based on the Goodenough-Kanamori 
rules of superexchange . ^^1^'' J. Goodenough in Ref. [H 
explained the FM sign of the nearest-neighbor coupling 
Ji in ^Cr2S(e)4 by indirect hopping between half-filled 
Cr t2g and empty eg states via p states of X anions. 
K. Dwight and N. Menyuk (Refs. [T9II20I ) analyzed vari- 
ous superexchange paths for interactions between up to 
6-th Cr neighbors and concluded that AFM coupling con- 
stants J3 between 3-rd neighbors and even weaker FM J4 
and Jg may be relevant alongside Ji . Then, the estimated 
J„ were used to examine the stability of different spiral 
ground states in ZnCr2Se4. To our knowledge the only 
attempt to obtain the values of J„ from ab initio band 
structure calculations was made in Ref. [2l| where the cou- 
pling constants J1-J3 were calculated by comparing the 
total energies of several simple spin configurations. 

The aim of the present work is to compare electronic 



2 



band structures of ACr2X4 {A=Zn, Cd, Hg and X=0, 
S, Se) spinels calculated within the local spin density 
(LSDA) as well as LSDA+C/ approximations and to es- 
timate exchange coupling constants between Cr spins by 
fitting the calculated energy of spin spirals to a classi- 
cal Heisenberg model. The paper is organized as follows. 
The spinel crystal structure is shortly described in Sec.lTTl 
Some details of the calculational procedure are given in 
Sec. mil In Sec. |TV] the band structures of the ACT2X4 
spinels calculated using LSDA and LSDA-|-[/ are com- 
pared and their dependence on the chemical composition 
is analyzed. The results on the exchange coupling con- 
stants Jn are presented in Sec. |Vl The comparison of 
calculated J„ to experimental data and the discussion of 
their origins are given in Sec. IVII Finally, the results are 
summarised in Sec. I VIII 



II. CRYSTAL STRUCTURE 

ACr2X4 compounds considered here belong to a large 
family of A^^ 32^^ Xl~ spinels which crystallise to a cubic 
Fd3m (N227) structure; with A, B, and X ions occupy- 
ing 8a (1/8,1/8,1/8), \M (1/2,1/2,1/2), and 32e (x,a;,x) 
Wyckoff positions, respectively. The experimental values 
of the lattice constant (oq) and fractional coordinates 
(x) are collected in Table HI The spinel crystal structure 
plotted in Fig. [T] can be considered as built of distorted 
Ct^X/^ cubes which share a Cr site. Cr and X ions that 
belong to the same cube form two regular tetrahedra with 
a common center, which coincides with the center of the 
cube. Each cube is linked via X ions to four regular 
AX/i tetrahedra. The centers of the Cr4X4 cubes and 
AX4 tetrahedra form two diamond lattices shifted by a 
vector (1/4,1/4,1/4). Finally, Cr ions are arranged along 
chains running in (110) directions and form the so-called 
pyrochlore lattice which consists of corner-sharing regu- 
lar tetrahedra. 

As the ionic radius of a cation increases in the 
row Zn^+ -^Cd^+ ^Hg'^^, X ions are pushed further 
away from the center of a AX/^ tetrahedron. This leads 
to lattice expansion and corresponding increase of Cr- 
Cr distances. The distortion of Cr4X4 cubes also in- 
creases, which allows to avoid too strong elongation of 
Ct-X bonds and manifests itself in the increase of the x 
parameter. When ions are replaced by S^~ or Se^~, 
which have significantly larger ionic radii, the length of 
both A-X and Cr-X bonds increases. This results in 
even stronger increase of the lattice constant. The near- 
est Cr-Cr (dcr-Cr), Cr-X {dcr-x), and A-X {(Ia-x) dis- 
tances are summarized in Table HI 

Each Cr site is surrounded by a trigonally distorted 
Xg octahedron; with all Cr-X distances in the octahe- 
dron being equal. The degree of the trigonal distortion 
is determined by the value of the fractional coordinate x. 
For X >0.25 the octahedron is expanded along one of the 
(111) directions and becomes regular for x=l/4. 

The local symmetry of a Cr site is D^d- It is worth 




FIG. 1: (Color online) The crystal structure of ACX2X4, 
spinels with an A ion shifted to the origin. Distorted Cr4X4 
cubes and AX4, tetrahedra are plotted by thick (red) and thin 
(green) lines, respectively. 



TABLE I: Experimental lattice constants ao (A) and frac- 
tional coordinate a; of X ions used in the calculations and 
the shortest Cr-Cr, Cr-X, and A-X distances (A). The last 
column contains experimental values of Ocw taken from Ref. 
i 



Compound 


a 


X 


dCr-Cr 




dA-X 


Ocw 


ZnCr204" 


8.327 


0.2616 


2.944 


1.990 


1.970 


-398 


CdCr2 04'' 


8.600 


0.2682 


3.041 


2.006 


2.133 


-71 


HgCraOi" 


8.661 


0.2706 


3.062 


2.003 


2.184 


-32 


ZnCr2S4" 


9.982 


0.2619 


3.529 


2.383 


2.367 


7.9 


CdCr2S4'' 


10.240 


0.2647 


3.620 


2.419 


2.478 


90 


HgCr2S4'= 


10.256 


0.267 


3.626 


2.402 


2.523 


140 


ZnCr2Se4^ 


10.484 


0.2599 


3.707 


2.522 


2.450 


155 


CdCr2Se4« 


10.735 


0.2642 


3.795 


2.540 


2.588 


184 


HgCr2Se4" 


10.737 


0.264 


3.796 


2.543 


2.585 


200 



"H. Sawada, Ref.IH 

''H. Ucda, Ref.|l0|. 

'^E. Riedel and E. Horvath, Ref. [H. 

''T. Borovskaya et al, Ref [pi. 

■^J. Hemberger et al., Ref. 

•'j. Akimitsu et al, Ref. Hi]. 

9J. Krok-Kowalski et al, Ref.IH. 

noting that since six Cr ions nearest to a Cr site form a 
trigonal antiprism, the symmetry remains trigonal even 
when X is equal to the ideal value of 1/4. 



III. COMPUTATIONAL DETAILS 

The calculations of the electronic band structure of the 
ACv2Xji spinels were performed for the experimentally 
observed lattice constant and X fractional coordinates 
(Table [l| using the hnear muffin-tin orbital (LMTO) 
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method^ with the combined correction terms taken into 
account. In order to decrease overlap between atomic 
spheres 3 sets of empty spheres (E) were added at 86 
(0,0,0), 16c (3/8,3/8,3/8), and 48/ (x', 1/8, 1/8) Wyckoff 
positions. Mufhn-tin orbitals with angular momentum 
I < 2 for A, Cr, and X spheres and I < 1 for empty 
spheres were included into the basis set. The Perdew- 
Wang parameterisation^ for the exchange-correlation 
potential in the local spin-density approximation was 
used. Brillouin zone (BZ) integrations were performed 
using the improved tetrahedron method.— 

Calculations for spiral spin structures based on the 
generalized Bloch theoreni^^ were performed under the 
assumption that the direction of the magnetisation m is 
constant within each atomic sphere. Inside a sphere at 
t + R, where t defines its position in a unit cell and R is a 
lattice vector, the magnetisation direction is determined 
by an angle 9t between the spin moment and z axis, the 
wave vector q of a spiral, and a phase c/jt'- 



following block form: 



m(rt) = m{rt) 



cos((/>t -I- q • R) sin( 
sin(^t -I- q • R) sin( 
cos 6t 



(1) 



with rt = r — t — R. 

When performing LMTO calculations for spin spirals 
it is convenient to split the LSDA exchange-correlation 
potential into a spin- and q-independent part V — {V^ + 
Vi)/2 and an effective exchange field B = {V^ — Vi)/2, 
where and Vi are exchange-correlation potentials for 
majority- and minority-spin electrons defined in the site- 
dependent local spin frame, in which the spin-density 
matrix for a given site is diagonal. In the present cal- 
culations spin-independent LMTO basis functions were 
constructed starting from the solution of the Kohn- 
Sham equation with only the spin-independent part V of 
the exchange-correlation potential included to the LSDA 
one electron potential. Matrix elements of the spin- 
dependent part of the LMTO Hamiltonian 



X 



t,R t 

COS St sin6lt^j^e-'(**+i-^) 



sin6it^P^e^('^'+i-^) 



' cos Ut 



,(2) 



where cr is a vector of Pauli matrices, were included at a 
variational step. Radial matrix elements of _B(r) between 
the solution of the Kohn-Sham equation inside a sphere 
(j)^{r) and/or its energy derivative tpuii") were calculated 
by numerical integration. 

The diagonal in spin indices matrix elements of Hb 
calculated between two Bloch wave functions with wave 
vectors k and k' do not depend on q and, as usual, are 
non-zero only if k = k'. In the absence of the spin-orbit 
coupling the only off-diagonal in spin indices terms of 
the Hamiltonian are Hg {q) and Hg\q) given by Eq. 

which couple the states with k' — k = ±q. Then, 
the LMTO Hamiltonian matrix can be written in the 



ii 



H 



-^k-q/2,k-q/2 



it 



k-q/2,k+q/2 
rrTT 

k+q/2,k-q/2 -"k+q/2,k+q/2 



(3) 



Thus, for a spin spiral with an arbitrary q one needs to 
diagonalise one 2N x 27V matrix instead of two N x N 
matrices, one for each spin channel, in conventional spin- 
polarised calculations.— If, however, the spin-orbit cou- 
pling term is included into the Hamiltonian it addition- 
ally couples majority- and minority-spin states with the 
same k and the Hamiltonian matrix becomes infinite. 

Most calculations reported in Sec. |V] were performed 
for planar spin spirals with all 9t — tt/2. Then, for a 
given q the magnetisation direction inside a sphere is 
defined solely by the phase (j)t. The phases at Cr sites 
were fixed by requiring that 0( = q • t. At other atomic 
and empty sites they were determined self-consistently. 
At each iteration of a self-consistency loop a rotation to 
the local spin frame in which the spin-density matrix be- 
comes diagonal was found and the corresponding rotation 
angles were used to determine new magnetisation direc- 
tion for the next iteration.— Iterations were repeated un- 
til self-consistency in the electron spin density as well as 
in the magnetisation direction in each atomic sphere was 
achieved. 

Finally, in order to account for correlation effects in the 
Cr 3d shell we adopted the LSDA-I-C/ method^^ in the 
rotationally invariant representation . ^^'^'^ The so-called 
atomic limit^ was used for the double counting term. 
The effective screened Coulomb repulsion U between 3d 
electrons was considered as a parameter of the model and 
varied from 2 to 4 eV. For the on-site exchange integral 
Jh the value of 0.9 eV estimated from the LSDA calcu- 
lations was used. 



IV. COMPARISON OF THE ELECTRONIC 
STRUCTURES OF 0-, S-, AND Se-BASED 
SPINELS 



Spin-restricted LDA results and hopping matrix 
elements 



The effect of the chemical composition on the elec- 
tronic structure of AGT2X4 spinels can be analyzed by 
comparing the densities of Cr d and X p states ob- 
tained from spin-restricted LDA calculations for ACr2X4 
(Fig. Since site-resolved densities of states (DOS) for 
the v4Cr2S4 spinels show the same trends as those for 
74Cr2Se4, only DOS calculated for ACr204 and ACr2Se4 
compounds are presented in Fig. [2] 

In all the ylCr2Ar4 spinels considered in the present 
work occupied bands in the energy range down to —10 
eV below the Fermi level {Ep) originate from A d, Cr d, 
and X p states. A d states in the sulfides and selenides 
as well as Cd 4d states in CdCr204 form a narrow group 
of bands below the bottom oi X p states. Zn 3c? and Hg 
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FIG. 2: (Color online) Densities of Cr d (red solid lines), X 
p (blue dashed lines), and A d (green dotted lines) states in 
ACT2O4 and j4Cr2Se4. Zero energy is chosen at the Fermi 
level. 



FIG. 3; (Color online) Symmetry-resolved densities of the Cr 
d aig (red solid lines), (blue dashed lines), and (green 
dotted lines) states in j4Cr204 and ylCr2Se4. Zero energy is 
chosen at the Fermi level. 



5d bands in the corresponding oxides cross the bottom of 
O p states and hybridise strongly with the latter which 
results in the appearance of wide peaks of d DOS between 
—9 and —7 eV. The A d states are completely occupied 
and lie well below Ep so that they have no effect on the 
magnetic properties of the ACT2X4 spinels. 

O p-derived bands in ACr204 spread over the energy 
range from —8 to —3 eV and are separated by a gap of 
^ 1.5 eV from Cr 3d states which give prevailing con- 
tribution to the bands crossing the Fermi level and a 
prominent DOS peak at Ep. Because of weaker elec- 
tronegativity of and Se^~ ions as compared to 0^~, 
X p states in ACr2S(e)4 spinels move closer to Ep and 
form bands between —6.5 and —0.5 eV. In contrast to 
oxides the top oi X p bands overlaps with the bottom 
of Cr d ones and the gap between X p and Cr d states 
closes. 

Cr d states are split by the cubic component of the 
crystal field at a Cr site into a triplet t2g {dxy, dzx, and 
dyz) and a doublet eg {dsz2__i and dx2-y2). In a CrXg 
octahedron the t2g and Cg states form relatively weak 
pdn- and much stronger pda-type bonds with the X p 
states, respectively. Bonding Cr d - X p combinations 
participate in the formation of the X p-derived bands 
which is evidenced by the rather high density of Cr d 
states in this energy range (see Fig. [5]). These states 
are completely filled and stabilise the GiXq octahedron. 
In all the ACr2X4 spinels considered here, the partially 
occupied bands crossing Ep are formed by antibonding 
combinations of the Cr d t2g and X p states with the 
dominant contribution of the former. These states, filled 



with 3 electrons, play the crucial role in the formation of 
Cr magnetic moments and effective exchange interactions 
between them. The Cr d states of the Cg symmetry which 
form antibonding pda combinations with the X p states 
are shifted to higher energies and separated by an energy 
gap from the antibonding Cr t2g - X p states. For the 
sake of brevity in the following these bands are referred 
to simply as Cr d t2g and eg bands. 

As the local symmetry of a Cr site is lower than cubic 
the t2g states are additionally split into a singlet and a 
doublet which transform according to aig and eg repre- 
sentations of the D^d group. We denote this doublet as 

in order to distinguish it from the doublet formed by 
d^z^^i and d^-i^y^ orbitals (e^) which transforms accord- 
ing to the same Cg representation. The densities of the 
Cr d states of aig, e^, and symmetries are shown in 
Fig. [3 

The trigonal splitting between the aig and states 
is much smaller than their bandwidths; with the cen- 
ter of gravity of the former being about 0.05 eV lower. 
The corresponding DOS curves in ^Cr204 are, however, 
remarkably different. The total width of the Cr t2g sub- 
band is determined by the states. Their DOS curve is 
very asymmetric with low DOS at the bottom and a huge 
DOS peak at the top of the t2g bands. The aig states 
are significantly narrower. They are responsible for the 
DOS peak at Ep but contribute also to the high-energy 
eg DOS peak. This strong suppression of the width of 
the aig DOS is a characteristic feature of 3d transition 
metal oxides with the spinel structure and a small ionic 
radius of an A ion. Recently, two scenarios based on a 
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strong Coulomb interaction within V d shell and the lo- 
calized behavior of V aig states have been proposed to 
explain heavy- fermion-like properties in LiV2 041^1^ 

A simplified analysis using Slater-Koster integrals^! 
shows that the narrowing of aig DOS is caused by a 
delicate balance between direct d-d hopping matrix el- 
ements t(id between the Cr t2g states and an effective in- 
direct hopping t^j^ via O p states. Under the assumption 
that tddS = tdda/G, tddrr = -2tdda/3, and tdda < the 
strongest hopping between t2g states of nearest Cr ions 
is a negative dda-type hopping between those orbitals 
which have their lobes directed along one of Cr chains, 
e.g. a dxy-dxy hopping along the [110] direction. Hop- 
pings of ddiT type between the other two t2g orbitals {dyz 
and dzx) along the same chain are positive and about 
two times weaker, whereas d^y-dyz/zx hoppings are zero. 
After making transformation to the symmetrised combi- 
nations of t2g orbitals it turns out that the matrix el- 
ements between aig and orbitals on neighboring Cr 
sites are much smaller than the matrix elements between 
the states of the same symmetry, with the aig~aig hop- 
ping being stronger than the Sg-^g one. This results in 
almost decoupled aig and bands, with the bandwidth 
of the former being larger than of the latter. 

The contribution of the hopping via O p states to the 
lowest order in f^^^/ied — Sp), where tpd is the Cr d-0 p 
matrix element and Sp and Sd are the energies of the Cr d 
and O p states, can be easily estimated using the Lowdin 
partitioning technique.'^- In the ideal {x = 1/4) spinel 
structure with regular Cr404 cubes the indirect hopping 
contributes only to dyz-dzx and dzx-dyz hybridisations 
via O Pz states. The corresponding matrix element t^z.xz 
is negative, so that the direct and indirect terms compen- 
sate each other. However, as it was shown for LiV204 in 
Ref. '39'and verified for ZnCr204 in the present work, the 
densities of aig and states calculated for a hypotheti- 
cal spinel structure with x — 1/4 are of the same width. 
Thus, the reduction of the dyz-dzx hopping is not suffi- 
cient to explain the narrowing of the aig states. When 
the distortion of Cr4X4 cubes is taken into account, non- 
zero indirect hoppings proportional to the degree of the 
distortion 6 ~ a; — 1/4 appear between all pairs of t2g 
orbitals. However, the most significant one is an indi- 
rect dxy-dxy matrix element between the orbitals point- 
ing along Cr chains. Already for x—0.26 it becomes as 
large as 0.6tPz.xz and, being of the opposite sign, strongly 
suppresses the direct dxy-dxy hopping. The net effect of 
the competition between the direct and indirect contri- 
butions to the effective d-d hopping is the reduction of 
the aig bandwidth relative to the width of states. If, 
however, the direct hopping becomes weak, the width of 
the aig and states is governed by the indirect hopping 
via X p states and depends strongly on the deviation of 
the fractional coordinate x from the ideal value of 1/4. 

Finishing the discussion of the d-d hybridisation we 
have to mention strong hopping matrix elements between 
t2g and Cg states which appear as a result of either direct 
hopping between, for instance, dxy and d^z^-i orbitals in 



the [110] direction or indirect one via Px and Py states of 
two X ions closest to the Cr-Cr bond. The indirect hop- 
ping is allowed already in an undistorted CT4X4 cube and 
is stronger than t^z xz for one of the involved p-d bonds 
is of pda type. However, the direct and indirect contri- 
butions are of opposite signs which may lead to strong 
suppression of effective t2g-eg hopping matrix elements. 
Finally, an indirect contribution to hopping between 
states of nearest Cr ions is proportional to 6 and is of the 
same sign as the corresponding direct hopping. 

Let us come back to the results of band structure cal- 
culations shown in Figs. and [3] As Zn is replaced by 
a larger Cd or Hg ion, the width of the Cr t2g states de- 
creases from 1.8 eV in ZnCr204 to 1.6 eV in HgCr204. 
The narrowing of the bands can be explained by the re- 
duction of the strength of both d-d and Cr d-0 p hop- 
ping matrix elements caused by the lattice expansion (see 
Table IJ). The shape of aig and DOS, however, also 
varies considerably. This indicates the change of relative 
strengths of the direct d-d and indirect d-0 p-d hoppings 
between Cr t2g states. 

Cr Cg states are affected much stronger by the change 
of the chemical composition because of their interaction 
with s bands of an A ion. In ZnCr204 bands cross just 
the bottom of a Zn 4s band near the BZ center. Cd 5s 
and, especially, Hg 6s bands shift to lower energies and 
overlap with Cr bands. In CdCr204 the hybridisation 
with a Cd 5s band is responsible for a low energy tail 
of Cr DOS. In HgCr204 the bottom of a Hg 6s band 
comes so close to Cr t2g bands that it starts to hybridise 
via O p states with Cr aig bands, which is evidenced by 
the appearance of noticeable density of aig and states 
in the energy gap between t2g and e^-derived bands. 

Symmetry resolved densities of Cr d states in selenides, 
presented on the right hand side of Fig. [31 are strikingly 
different from the corresponding DOS curves in oxides. 
The densities of aig and states do not differ so much 
as in oxides; with the top of the bands being shifted 
to somewhat higher energies. The huge DOS peak found 
at the very top of the t2g states in ZnCr204 and slightly 
lower in the two other oxides shifts in selenides closer 
to the center of the t2g bands. The energy difference 
between the Cr and t2g states is somewhat smaller 
than in oxides. The DOS curves become wider and 
more symmetric. Another important distinction of the 
ACr2S(e)4 electronic structure is the enhanced weight of 
the states in the Bloch wave functions of Cr t2g bands, 
which is revealed by noticeable density of the states 
in the corresponding energy range. 

The difference between the dispersions of Cr d bands 
in the oxides and selenides is illustrated in Fig. 0] which 
shows "fat" bands calculated for ZnCr204 and ZnCr2Se4; 
with the size of filled circles and squares being propor- 
tional to the partial weight of Cr aig and states, re- 
spectively, in the Bloch wave function at a given k-point. 

The abovementioned differences in the band disper- 
sions and DOS suggest that the effective d-d hopping 
matrix elements should also be significantly different. 
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FIG. 4: (Color online) "Fat" bands calculated along some 
symmetry lines in the fee BZ for ZnCr2 04 (lower panel) and 
ZnCr2Se4 (upper panel). The size of filled (red) circles and 
(blue) squares is proportional to the partial weight of Cr aig 
and eg states in the Bloch wave function, respectively. Zero 
energy is chosen at the Fermi level. 



which is not surprising taking into account the increase 
of the lattice constants in ACr2S(e)4 as compared to 
ACr204 (Table [H. The increase of Cr-Cr distances leads 
to strong reduction of the direct d-d hopping matrix el- 
ements which fall off as l/dQj._Q^.. It should be noted, 
however, that a tight binding analysis of the band struc- 
ture of CdCr2S4 performed in Ref. I40I showed that the 
direct hopping between Cr d states is not negligible even 
in sulfides and selenides. The contribution of the indi- 
rect hopping via X p states to the effective d-d hoppings, 
which is proportional to t^^Ksd — £p), does not change 
much because i) S ip and Se 4p states are significantly 
more spatially extended than O 2p states and, although 
dcr-x also increases in y4Cr2S(e)4, Cr d~X p hoppings tpd 
remain approximately of the same strength as in ylCr2 04; 
ii) the energy difference £4 — Sp in sulfides and selenides 
is smaller than in oxides. 

Similar to oxides the Cr t2g states in selenides become 
narrower with the increase of the A ionic radius. The 
hybridisation of Cr states with Hg 6s bands leads to 
sharpening of a DOS peak near the bottom of the 
sub-band. 




-8-6^-2 2 4 -8 -6 -4 -2 2 4 



Energy (eV) Energy (eV) 

FIG. 5: (Color online) Cr d (red solid lines) and X p (blue 
dashed lines) DOS in ylCr204 and ACr2Se4 obtained from 
spin-polarized calculations with FM alignment of Cr magnetic 
moments. Zero energy is at the Fermi level. 



B. The effect of spin-polarisation 

When in spin-polarized LSDA calculations Cr d states 
with two spin projections are allowed to have different oc- 
cupations, the strong on-site exchange interaction splits 
the half-filled Cr t2g shell into occupied majority-spin i2gT 
and unoccupied minority-spin t2gx- The densities of Cr 
d and X p states in ^Cr204 and ACr2Se4 obtained from 
spin-polarized calculations with the ferromagnetic (FM) 
alignment of Cr moments are shown in Fig. [S] 

In all the spinels studied here the exchange splitting 
of Cr d states is about 2.7 eV which gives an estimate of 
0.9 eV for the on-site Hund's exchange coupling Jh- The 
exchange and crystal field splittings are of comparable 
strengths, so that Cr t2gi states are found at the same 
energy as e^^ . In the oxides t2g^ states remain separated 
by an energy gap from O p bands. In ZnCr204 and 
CdCr204 the FM solution is insulating with completely 
occupied t2gi and empty t2gi states. In IIgCr204 the 
bottom of a minority-spin Hg s band hybridised with 
Cr aigi states crosses t2g] bands leading to a metallic 
solution. 

In contrast to the oxides, in the ^Cr2S(e)4 spinels t2gi 
states move below the top of the X p states. As a re- 
sult the highest occupied majority-spin bands are formed 
mainly hy X p states strongly hybridised with Cr e^^ . In 
the Zn and Hg spinels they overlap with the very bottom 
of mostly unoccupied i2gx bands while for CdCr2S4 and 
CdCr2Se4 insulating solutions with a tiny gap were ob- 
tained. The LSDA band structures and DOS calculated 
for CdCr2S4 and CdCr2Se4 are in good agreement with 
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the results of previous calculations j i In the present 
work, however, the FM solution for CdCr2S4 is insulating 
whereas in Ref. 113 metallic FM solutions were obtained 
for both compounds. A possible cause for the difference 
is that in Ref. 40^ structural data with the S fractional 
coordinate x=l/4 were used. 

A Cr spin magnetic moment defined as a volume in- 
tegral of the spin density over a Cr atomic sphere varies 
from 2.86 /^b in ZnCr204 to 3.15 /ib in CdCr2Se4, which 
agrees with a good accuracy with the value of 3 fiB ex- 
pected for the completely polarized t2g shell filled with 
3 electrons. The hybridisation with Cr d states induces 
weak negative spin-polarisation oi X p states. It is, how- 
ever, partially compensated by a positive contribution 
coming from the "tails" of Cr d states that have the 
d symmetry inside an X atomic sphere. The net spin 
magnetic moment Mx induced on X ions is small and 
negative with Mq « —0.01 //b, Ms « —0.05 ^b, and 
Msc ~ —0.09 /iB. The moment induced on A ions is less 
than 0.05 /ib. 

The self-consistent FM solution discussed above can be 
obtained for any of the 9 spinels but it does not necessar- 
ily have the lowest total energy. If the coupling between 
Cr spins is antiferromagnetic (AFM) then, because of ge- 
ometrical frustrations intrinsic to the pyrochlore lattice 
formed by Cr ions, any configuration of classical spins 
such that the sum of 4 spin vectors sitting in the cor- 
ners of each tetrahedron is equal to zero is the ground 
state of the Heisenberg Hamiltonian with nearest neigh- 
bor interactionsi^ The requirement of zero total spin in 
each tetrahedron is satisfied, for example, if Cr spins are 
aligned ferromagnetically along the [1±10] chains and an- 
tiferromagnetically along [01±1] and [10±1] ones or, in 
other words, if Cr spins in each (001) plane are parallel 
but the spins in consecutive planes along the [001] direc- 
tion are antiparallel to each other. In the following this 
spin configuration will be denoted as AFZ. It is easy to 
check that in this case only 4 of 6 Cr-Cr bonds in each 
tetrahedron are AFM while the other two remain FM 
and the magnetic energy cannot be minimised for all 6 
bonds simultaneously if the coupling between Cr spins is 
antiferromagnetic. 

Another spin configuration which gives zero total spin 
in each tetrahedron and does not break the cubic sym- 
metry of the lattice is a non-coUinear one with Cr spins 
directed along the lines passing through the center of a 
tetrahedron, i.e., along one of the (111) directions; with 
all four Cr spins pointing either to or away from the 
center. Self-consistent band structure calculations per- 
formed with this configuration of Cr spins gave for all the 
ACT2Xi spinels insulating solutions with zero net mag- 
netic moment, which in the following is referred to as 
a ZM solution. Cr magnetic moments calculated in the 
local spin frame are about 0.1 /ib smaller than for the 
corresponding FM solution but still close to 3 /j,b. De- 
pending on the magnetisation directions of four Cr ions 
in a 014X4^ cube the magnetisations in four X atomic 
spheres are parallel or antiparallel to the vectors point- 



ing to the center of the cube so that small X magnetic 
moments also cancel each other. With this arrangement 
of Cr and X moments the spin moments of A ions are 
equal to zero. It is worth noting that for all the AGT2X4 
compounds the LSDA total energy difference between the 
AFZ and ZM solutions does not exceed 3 meV per for- 
mula unit. 

The comparison of the LSDA total energies of FM and 
ZM solutions shows that in ZnCr204 the solution with 
zero net moment is significantly more favourable; with 
the energy difference per formula unit (f.u.) being 171 
meV. In CdCr204 this difference decreases to 19 meV, 
whereas in HgCr204 the two solutions are almost degen- 
erate; the FM one being 2 meV lower. In the sulfides and 
selenides the FM solution becomes more preferable. As 
Zn is replaced by a heavier ion the energy difference be- 
tween the two solution increases from 43 to 118 meV/f.u. 
in sulfides and from 80 to 134 meV/f.u. in selenides. 

These results indicate that the change of the sign of 
dominant exchange interactions from AFM in oxides to 
FM in sulfides and selenides is captured already by LSDA 
band structure calculations. Moreover, all the com- 
pounds show a clear tendency to ferromagnetism as the 
lattice expands due to the increase of the radius of A 
ions. 



C. LSDA+C/ results 

In ACT2X4, as well as in many other 3d compounds, 
LSDA underestimates correlation effects in the rather lo- 
calised Cr 3d shell. The strong electronic correlations can 
be accounted for, at least on a mean-field level, by using 
the LSDA-|-[/ approach. Since Cr t2g states are split by 
the on-site exchange interaction into occupied majority- 
and unoccupied minority-spin states, charge and orbital 
degrees of freedom are frozen already in LSDA. As a 
consequence, LSDA-|-f7 band structures of the 74Cr2Ar4 
spinels are qualitatively similar to the LSDA results. 

When non-spherical terms of the Coulomb repulsion 
are neglected, the orbital dependent LSDA+U potential 
Vi can be approximated as Vi = C/'(l/2 — n^), where Ui is 
the occupation of i-th localised orbital, U' — U~Jyi, and 
U is the effective screened Coulomb repulsion between 3d 
electrons. Taking into account that the orbital occupa- 
tions rii of Cr t2gT ^^'^ ^2gi states are close to 1 and 0, 
respectively, the main effect of the LSDA-I-J7 potential is 
to shift the t2g| states to lower and t2gi to higher ener- 
gies by U' /2. Thus, the splitting between the t2gi and 
t2gi states increases from 3Jh in LSDA calculations to 
3Jh + U' in LSDA+C/ ones. In oxides, because of the 
downward shift of the t2g| states, they start to overlap 
with O p bands. In v4Cr2S(e)4 the t2g^ states move fur- 
ther below the top of X p bands and their contribution to 
the highest occupied majority spin bands decreases. Due 
to the increase of the t2gX~i2gT splitting, ferromagnetic 
LSDA-f[/ solutions for all 9 compounds considered in 
the present work become insulating starting from U =3 



8 



eV. 

Because of the strong Cr e'^-X p hybridisation dis- 
cussed in Sec. lIVAl the orbital occupations rii calculated 
for the formally unoccupied Cr states are in the range 
0.28-0.48 for the majority- and about 0.22-0.25 for the 
minority-spin orbitals. As a result the shift of the Cr 
states to higher energies caused by the Coulomb repul- 
sion U is somewhat smaller than for the t2gi states for 
which Hi does not exceed 0.07 and the splitting between 
and t2gi states decreases with the increase of U. 



V. EFFECTIVE EXCHANGE COUPLING 
CONSTANTS 

Effective exchange coupling constants J„ between Cr 
spins were determined by performing band structure cal- 
culations for spiral spin structures with a varying wave 
vector q. Then, the q-dependence of the energy of 
the spin spirals was mapped onto a classical Heisenberg 
model. In order to get reliable estimates for J„ one needs 
to perform calculations for a sufhciently large number of 
q. Since self-consistent calculations for spin spirals are 
much more time-consuming than conventional coUinear 
spin-polarised calculations, the q-dependence of their en- 
ergy was calculated using the so-called local force the- 
orem (LFT)42ii^ According to this approach the total 
energy difference between two spin configurations can 
be approximated by the difference of their band ener- 
gies, provided that the calculations are performed start- 
ing from the same electron spin densities and only the 
magnetisation direction varies. Strictly speaking this ap- 
proximation is justified only for small deviations of the 
magnetisation direction from some collinear spin arrange- 
ment, i.e., for small |q|, and before using it for short wave- 
lengths we have numerically checked its accuracy for the 
case of the ACr2X4 compounds. 

First, self-consistent calculations were performed for 
two sets of planar spin spirals with the magnetisation di- 
rection inside a Cr sphere fixed by polar angles 6 — n/2 
and = q • (t -I- R) , where t is the position of a Cr site in 
the unit cell and R is a lattice vector. The wave vectors 
q=(0,0,g) and q=(g,q,0) of the spirals varied along the 
T-X and T~K high symmetry directions, respectively, in 
the range < g < 2 in 27r/ao units. Obviously, the (0,0,0) 
spiral is the collinear FM structure with the magnetisa- 
tion directed along the x axis, whereas (2,2,0) and (0,0,2) 
spirals are equivalent to the collinear AFZ spin structure 
with zero net magnetisation discussed in Sec. IIVBI The 
magnetisation directions in other atomic spheres were de- 
termined self-consistently. 

Then, band energies for the same spin spirals were cal- 
culated starting from self-consistent FM electron densi- 
ties for ACr2S(e)4 and from the densities for the ZM so- 
lution in the case of oxides. When using LFT for spinels 
one needs to specify the magnetisation direction not only 
for magnetic Cr ions but also in A, X, and E spheres in 
which small but finite magnetic moments are induced. 
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FIG. 6: (Color online) The dependence of the total Etoti'i) 
and band Ebnd{<\) energies calculated within the LSDA on the 
wave vector of (0,0, g) and (q,q.ff) spin spirals. The energies 
are plotted relative to the energy of the FM states with g=0. 



In a general case this direction is not uniquely deter- 
mined by q but instead depends in a nontrivial way on 
the orientation of Cr magnetic moments. We found that 
good agreement between the q-dependencies of the total 
-Etot(q) and band £'b„d(q) energies can be obtained when 
the small exchange splitting vn A, X, and E spheres is 
completely neglected. Ebndil) curves calculated in this 
way for some of the Cr spinels are compared to Etot (q) h^ 
Fig. [6l Both energies are plotted relative to the energy 
of the FM solution with q=0. The difference between 
Ebndi<l) and Etot{(l) does not exceed 5 meV and is much 
smaller than the variation of the total energy with q even 
for CdCr204. Somewhat larger differences of about 8 
meV are calculated for sulfides (not shown in Fig. 
If the spin-polarisation inside non-Cr spheres was not 
switched-off and the magnetisation directions obtained 
from self-consistent calculations were used in LFT calcu- 
lations, slightly larger differences between the total and 
band energies were obtained. Below, when discussing ex- 
change coupling constants between Cr spins, we will drop 
the subscript bnd and denote the band energy calculated 
with the spin-polarisation in A, X, and E spheres set to 
zero simply as -E'(q). 

The calculations for spin spirals confirm the conclu- 
sion, drawn in Sec. IIVBI on the base of comparison of 
the FM and ZM total energies, that the energy of the 
zero net moment solution {q=2) is lower than the FM 
one (<z=0) in oxides whereas in ACr2S(e)4 the FM solu- 
tion becomes more stable. However, the shift of the min- 
imum of E{q) curves to intermediate values of q, most 
clearly seen for the T-K direction, witnesses that along- 
side geometrical frustrations there is also a competition 
between the nearest-neighbor exchange interaction and 
more distant ones. Only in CdCr2Se4 and HgCr2Se4 the 
minimum is at g=0 and the FM solution has the lowest 
energy. 

Quantitative estimates for J„ up to the 4-th shell of Cr 
neighbors were obtained by mapping the energy of spin 
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FIG. 7: (Color online) Exchange coupling constants J„ up to 
the 4-th Cr neighbors. 

spirals onto a classical Heisenberg model 

H 4 4 z„ 

^H^^EEE-^-s.-s,, (4) 

i—l n—1 j—1 

where i numbers Cr sites in the unit cell and j runs over 
Zn neighbors in the n-th shell around the site i. Positive 
J„ correspond to AFM coupling between Cr spins. Note 
that an additional factor of 1/2 appears in Eq. ^ be- 
cause the magnetic energy is given per formula unit and 
there are two formula units in the unit cell. The length 
of the Cr spins was fixed to 3/2. Then, -EH(q) de- 
pends only on the angle between a pair of Cr spins which 
is uniquely determined by the wave vector q of a spiral 
provided that for 4 Cr ions at positions in the unit cell 
the phases are fixed by 0i = • q. 

Some pairs of Cr ions coupled by J1-J4 are shown in 
Fig. [71 Each Cr site has 6 first (Ji) and 12 second ( J2) 
neighbors at distances aoV2/4 and ao\/6/4, respectively. 
The 3-rd shell consists of 12 Cr sites at ooa/2/2 which 
are split into two inequivalent sextets. The correspond- 
ing coupling constants are denoted as J3 and J3 . J3 
couples Cr sites which lie on one of (110) chains and are 
actually the 2-nd Cr neighbors along the chain. Cr sites 
coupled by J3 belong to parallel Cr chains. For the spin 
spirals considered here terms proportional to J3 and J3 
in Eq. ^ have the same q-dependence so that it was not 
possible to separate their contributions to the magnetic 
energy and only their average J3 = (Jg -I- J3 )/2 could be 
determined from the fit. We will return to the discussion 
of J3 and J3 later. Finally, in the 4-th shell there are 
12 Cr sites at the distance ao\/T0/4 (J4) which lie at the 
same chains as the 2-nd Cr neighbors. 

In order to get reliable values of J„, additional LFT 
calculations for spirals with wave vectors {q,q,q), (1,9,0), 



(0,0,t« (q.q.O) 




q (2ii/a) q (2ii/a) 

FIG. 8: (Color online) The comparison of calculated band 
energies (i?bnd) and the results of the fit using Eq. Q (-Eflt) 
for q=(0,0,g) (left panels) and q=(g,g,0) (right panels). Indi- 
vidual contributions En of the terms proportional to Jn are 
also plotted. 



and (1,1, <;) were performed. For the (l,g,0) and (1,1, q) 
spirals the contributions to E^h (q) proportional to J3 and 
Ji, respectively, do not depend on q, which allows to 
determine these coupling constants with higher accuracy. 
Explicit expressions for -Eh (9) for the above mentioned 
spin spirals are given in Appendix. 

The results of a simultaneous least-squares fit of Eb_{ci) 
given by Eq. (|4]) to E{q) calculated within LSDA along 
the five q directions are shown in Fig. [H] together with 
contributions to Eii{q) coming from different nearest 
neighbors shells (En). The values of J„ obtained from 
the fit are collected in Table [III The dominant con- 
tribution to £'H(q) is provided by the nearest neigh- 
bor coupling Ji which is AFM in ACr204 and FM in 
ylCr2S(e)4. In oxides, however, the strength of AFM Ji 
decreases dramatically with the increase of the A ionic 
radius: in CdCr204 Ji is more than 5 times weaker than 
in ZnCr204, while in IIgCr204 it becomes almost zero. 
In ^Cr2S(e)4 the strength of FM Ji tends to increase in 
the row Zn^Cd^Hg but the changes are not as strong 
as in the oxides; with the values of Ji calculated for the 
Cd and Hg compounds being comparable. 

Another significant contribution to En^q) comes from 
the term proportional to J3 which is AFM in all the 
ACT2Xi spinels considered in the present work. It is the 
competition between the Ji and J3 terms which shifts 
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the minimum of the (0,0, q) and (q,qfl) curves to incom- 
mensurate q vectors. For q—{0,Q,q) Eii{q) has the same 
g-dependence as a Hnear chain of classical spins with 
competing nearest J ~ Ji and next-nearest J' = 2J3 
neighbor interactions. The ground state of such a chain 
becomes incommensurate if the ratio J'/ J is larger than 
the critical value of 0.25. Since J3 is not very sensitive to 
the size of A ion, the ratio 2J3/J1, which is less than 0.1 
in ZnCr204, increases to 0.45 in CdCr204, while in the 
Hg oxide J3 becomes the dominant magnetic interaction. 
The values of J3 calulated for 74Cr2S(e)4 are more than 2 
times larger than in the oxides. The largest \2J^/ Ji\ ra- 
tios of 0.75 and 0.48 are found for ZnCr2S4 and ZnCr2Se4, 
in which the FM nearest neighbor interaction is the weak- 
est. In other two selenides \2J^/ Ji \ « 0.26 is only slighly 
larger than the critical value of 0.25. 

According to the results of the the least-square fit 
the exchange interactions between the 2-nd and 4-th Cr 
neighbors are much weaker than Ji and J3. While the 
sign of J2 varies form one compound to another, J4 is al- 
ways ferromagnetic and somewhat larger in sulfides and 
selenides compared to oxides. 

The reliability of the fitted J„ was checked by compar- 
ing the values obtained from the fit to the band energies 
calculated for all 5 spin spirals {Jn^) to those fitted to 
(0,0,(7) and (g,<z,0) results only (J„ '). In the case of Ji, 
the relative uncertainty of the determination of J„, de- 
fined as 5n — \Jn^ — Jii^l Jn \ is Icss than 1%, except for 
HgCr204 where Ji itself is very small. The uncertainty 
in J3 values is about 1% in ACr2S(e)4 and 10% in oxides. 
Tha values of 82 ~ 5% and 84, ~ 10% are also smaller in 
ACr2S(e)4 than in the oxides, in which they are about 
30%. The largest uncertainties in the LSD A values of 
J2-J4 were found in the case of ZnCr2 04. A plausible 
reason is that the nearest-neighbors terms of the order 
of (Si • Sj)^ in the expansion of the magnetic energy^ 
may become comparable to the contribution of weak cou- 
plings between distant Cr neighbors. Neither the fitted 
values of J„ nor the quality of the fit were noticeably af- 
fected if exchange couplings between 5-th Cr neighbors 
were included into the fit. Accounting for more distant 
6-th neighbors allowed to decrease the mean deviation 
between the calculated and fitted energies for ACr2S(e)4 
but the obtained values of Jg were even smaller than J2 
and J4 and could not be determined reliably. Finally, sig- 
nificantly smaller values of (5„ were obtained when fitting 
i?H(q) to LSDA-I-C/ band energies. 

Jn obtained by fitting £'H(q) to band energies calcu- 
lated using the LSDA-I-C/ approach with U=2, 3, and 4 
eV are presented in Table IIIII Both Ji and J3 derived 
from the LSDA-I-C/ calculations with C/=2 eV for ACr204 
are weaker than the corresponding LSDA values. The 
only exception is IIgCr204, in which Ji is weakly FM 
in LSDA and becomes somewhat stronger in LSDA-I-C/. 
In ZnCr204 the value of Ji rapidly decreases with the 
increase of U as it is expected if J ~ /U . In Cd and 
Hg spinels, on the other hand, the increase of U tends 
to make Ji more ferromagnetic. It looks as if there were 



TABLE II: Exchange coupling constants Jn/fce (K) and Qcw 
(K) obtained from the least-squares fit of -E'H(q) given by Eq. 
Q to the energy of spin spirals calculated within LSDA. 





Ji/kB 


^2/^3 




J^/k-B 


6cw 




1 flQ 


1 s 




n 1 


Q1 (S 


CdCr2 04 


18 


-0.3 


3.8 


-0.9 


-172 


HgCr2 04 


-1 


1.8 


5.4 


-0.8 


-92 


ZnCr2S4 


-37 


2.8 


13.8 


-1.5 


49 


CdCr2S4 


-74 


0.6 


12.0 


-1.7 


392 


HgCr2S4 


-86 


2.8 


13.3 


-1.7 


432 


ZnCr2Se4 


-54 


0.9 


13.0 


-2.2 


228 


CdCr2Se4 


-88 


-0.7 


11.7 


-2.1 


526 


HgCr2Se4 


-86 


0.0 


11.6 


-2.6 


509 



two competing contributions of the opposite signs to Ji: 
as the AFM one is suppressed by U the ferromagnetic 
contribution wins. Comparable values of Ji=0.5 meV, 
J2 ~0 meV, and J3 = 0.15 mcV were obtained from 
LSDA-hC/ calculations for CdCr2 04 in Ref. |2l|. These 
calculations, however, were performed using theoretical 
lattice parameters. 

In ZnCr2S4 |Ji| calculated with C/=2 eV is slightly 
larger than the LSDA value and continue to increase with 
the increase of U , whereas in ZnCr2Se4 Ji is practically 
independent of U . In other ACr2S(e)4 spinels switching 
on U suppresses Ji but the increase of U from 2 to 4 eV 
has only minor effect on its strength. 

In all the compounds J3 obtained from the LSDA-I-C/ 
calculations is AFM and weaker than in LSDA. In 
ylCr204 and ^Cr2S4 its value decreases with the increase 
of U while in selenides its U dependence is very weak. J2 
remains vanishingly small but shows clear tendency to 
become more FM as U increases. Finally, in ACi20i and 
v4Cr2Se4 the increase of U affects J4 in opposite ways: in 
oxides the strength of FM J4 increases whereas in se- 
lenides I J4I decreases and for C/=4 eV it changes sign. 

In order to separate J3 and J3 contributions to En (q) 
we performed calculations for non-coUinear spin super- 
structures in which Cr moments only in every second 
(001) plane, the one that contains [110] Cr chains, are 
oriented in the ah plane; with their directions being de- 
termined as before by 9 — 7r/2 and ^ = q • (t -I- R). Cr 
moments in other planes, which contain [110] chains, are 
aligned parallel to the c axis {9 = 0) so that their direc- 
tions do not depend on (p and, consequently, on q. In 
this case the 3-rd Cr neighbors in the ab plane that lie on 
the same [110] chain ( J3) and on parallel chains ( J3') give 
different contributions to the g-dependence of the energy 
of {q,q,0) and {q,~qfl) spin spirals, which allows to de- 
termine J3 and Jg separately provided that other J„ are 
known. Calculations performed for ylCr204 and ACr2Se4 
show that Jg and Jg are of comparable strengths. The 
calculated J3 /J3 ratio is 0.6, 0.5, and 1.1 in Zn, Cd, and 
Hg oxides, respectively, while for the selenides the val- 
ues of 0.8, 0.8, and 1.0 were obtained. It seems plausible 
that the contribution of Hg 6s states to Cr d~X p~Cr d 
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TABLE III: Exchange coupling constants Jn/k^ (K) and 
Ocw (K) obtained from the least-squares fit of -EH(q) given 
by Eq. © to the energy of spin spirals calculated using the 
LSDA+t/ approach with U=2, 3, and 4 eV. 



ACr2X4 


u 








,/i / h-D 




ZnCr^O/i 


2 eV 


61 


0.3 


2.9 


-0.3 


-500 




3 eV 


40 


0.0 


2.3 


-0.4 


-328 




4 eV 


25 


-0.2 


1.9 


-0.5 


-209 


CdCroOi 


2 eV 


6 


-0.3 


2.2 


-0.4 


-64 




3 eV 


-4 


-0.4 


1.7 


-0.4 


12 




4 eV 


-9 


-0.4 


1.4 


-0.4 


62 




2 eV 


-7 


1.2 


3.3 


0.0 


-14 




3 eV 


-14 


0.6 


2.4 


-0.1 


59 




4 eV 


-18 


0.3 


1.9 


-0.2 


104 


/in ((To S/i 


2 eV 


-43 


1.7 


8.5 


-0.6 


175 




3 eV 


-48 


0.8 


7.2 


-0.8 


267 




4 eV 


-52 


0.1 


6.4 


-0.7 


306 


CdCroSi 

V^' ± ^^41: 


2 eV 


-62 


0.4 


7.2 


-0.9 


367 




3 eV 


-65 


-0.2 


6.2 


-0.9 


416 




4 eV 


-66 


-0.6 


5.7 


-0.7 


433 




2 eV 


-72 


1.6 


8.0 


-0.9 


412 




3 eV 


-74 


0.5 


6.6 


-0.9 


464 




4 eV 


-74 


-0.2 


5.8 


-0.7 


481 


ZnCr2Se4 


2 eV 


-49 


0.5 


8.3 


-0.9 


246 




3 eV 


-52 


-0.5 


7.8 


-0.6 


286 




4 eV 


-52 


-1.1 


8.1 


0.0 


286 


CdCr2Se4 


2 eV 


-70 


-0.5 


7.1 


-0.9 


432 




3 eV 


-69 


-1.1 


7.1 


-0.4 


434 




4 eV 


-67 


-1.4 


7.6 


0.2 


410 


HgCr2Se4 


2 eV 


-69 


-0.1 


7.5 


-0.9 


422 




3 eV 


-68 


-1.0 


7.2 


-0.2 


420 




4 eV 


-66 


-1.4 


7.8 


1.1 


384 



hybridisation is responsible for somewhat larger values of 
Jg in the Hg compounds. 

It should be mentioned that the values of J„ presented 
in Table Hill are somewhat different from the preliminary 
results published in Ref. '45' for the following reasons: i) 
for CdCr2 04 old structural data from Ref. with the 
O fractional coordinate a;=0.260, which seems to be too 
small, were used; ii) the calculations for oxides in Ref. 
I45I were performed with O s and p states only included 
in the LMTO basis set. This gives smaller values of Ji 
compared to the present calculations in which O d are 
also included into the basis; iii) the calculations in Ref. 
I45I were performed for [q^qfi) and (0,0,(7) spirals only and 
the coupling constants between 3-rd Cr neighbors along 
the (110) chains were included into the fit instead of J4. 
Later it was verified that accounting for true J4 to 4-th 
Cr neighbors improves the quality of the fit and seems 
more physical. We have to stress, however, that despite 
the differences in calculated values of J„ the main con- 
clusions of Ref. [4^ concerning the relative strengths of 
various exchange interactions and their origins remain 
unaffected. 



VI. DISCUSSION 

The effective Curie- Weiss temperatures estimated from 
the calculatated J„ according to 

C>CW = g 2^ ZnJn , (5) 

n 

with S=3/2 are given in the last column of Tables |TT] 
and mil The comparison of the estimated Ocw to ex- 
perimental values from Ref. i2 (Table |T| shows that the 
LSDA calculations reproduce the opposite signs of 9cw 
in ACr2 04 and ACr2S(e)4 but strongly overestimate its 
magnitude, which is not surprising as one can expect that 
LSDA underestimates correlation effects in the rather lo- 
calised Cr 3d shell. Since the absolute values of 8cw 
calculated for oxides rapidly decrease with the increase 
of U in the LSDA+C/ calculations, it is possible to adjust 
the value of U so as to reproduce the experimental Qqw- 
This, however, would require to use different values of 
U for the Zn, Cd, and Hg oxospinels. In ACr2S(e)4 the 
estimated values of 6cw: which are too high already in 
LSDA, further increase when the LSDA-I-J7 approach is 
used. As the dominant contribution to &cw is provided 
by Ji, it seems that the tendency to ferromagnetic cou- 
pling between nearest Cr neighbors is overestimated in 
the LSDA-I-J7 calculations. 

The comparison of the calculated J„ to experimental 
estimates for the exchange coupling constants in ZnCr2S4 
(Ji=2.66 K, J2=-1.15 K, J3=0.29 K, J4=0.13 K>12 
and CdCr2S4 (Ji=13.25 K, J3=-0.915 K)^ also shows 
that the calculations overestimate the strength of the ex- 
change interactions between Cr spins. Note that the ex- 
perimental values of Jn should be multiplied by a factor 
of 2 because of the different definition of the magnetic 
energy [see Eq. (2) in Ref. Il3. w hile theirs signs should 
be inverted since in Refs. [l2l47l the positive sign of J„ 
corresponds to FM coupling. Also it is worth mention- 
ing that Jn are not directly measurable quantities and 
in order to determine them from experimental data ad- 
ditional assumption should be made. The experimental 
wave vector qexp =(0,0,0.79) of the helimagnetic struc- 
ture in ZnCr2S4 is very close to qmm =(0,0,0.82) at which 
the LSDA energy of the (0,0,(7) spin spiral has a mini- 
mum. Under assumption that J2 and J4 are much smaller 
than Ji and J3 this gives the J1/J3 ratio of 0.39 and 0.37 
in the experiment and theory, respectively. For ZnCr2Se4 
LSDA gives the minimum of E{q) along the T-X at 
Qinm =(0,0,0.65) which is larger than the experimental 
wave vector q^xp =(0,0,0.466)." In the LSDA+C/ calcu- 
lations, however, the minimum shifts towards smaller q 
values. 

In order to understand better the dependence of the 
calculated exchange coupling constants on the chemi- 
cal composition and on the value of the parameter U 
in LSDA+C/ calculations, let us for simplicity consider a 
pair of Cr'^"'" ions with occupied t2g| states at zero en- 
ergy, unoccupied t2gi states at the energy A||, and an 
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FIG. 9: (Color online) A sketch of bonding and antibonding 
states between Cr tig state for FM (left) and AFM (right) 
alignments of Cr magnetisations. The Fermi level is plotted 
by a horizontal dotted line. 



FIG. 10: (Color online) A sketch of bonding and antibonding 
states between an occupied Cr t2g\ state and an unoccupied 
fig state with the same spin projection for FM (left) and AFM 
(right) alignments of Cr magnetisations. The Fermi level is 
plotted by a horizontal dotted line. 



effective tig-tig hopping matrix element t. The split- 
ting A|| is close to 3Jh in the LSDA and 3 Jh + U' in 
the LSDA+[/ calculations. Cr states with both spin 
projections are unoccupied and also split by A'^ which 
in LSDA is equal to 3Jh. The LSDA energy difference 
Ae_t between and tig states with the same spin pro- 
jections is equal to the crystal field splitting IS.cf which 
is of the same order as A'^^ as is evidenced by the results 
presented in Sec. IIVBI 

Hybridisation, either direct or indirect via X p states, 
between the half-filled Cr tig states results in AFM ex- 
change coupling provided that the tig bandwidth is less 
than A||. Indeed, if Cr moments are aligned ferromag- 
netically, the hybridisation between the tig states gives 
no gain in the kinetic energy as in the majority-spin 
chanel both bonding and antibonding combinations with 
energies = —t and e"'' = t are occupied, whereas in 

the minority-spin chanel both of them (e|'"'' — A|| ± t) 
are empty (see Fig. If the magnetisations of two Cr 
ions are antiparallel, an occupied tig^ state can hybridise 
only with the unoccupied tig state with the same spin 
projection of another ion lying at the energy A||. In 
this case only their bonding combination with the en- 
ergy £| « — i^/A|| is occupied. Since the energy of the 
bonding state with the opposite spin projection is exactly 
the same, this gives 

EpM — Eafm = '^JafmS'^ ~ 2i^/ A-^ (6) 

for the energy difference of the FM and AFM configura- 
tions. 

Cr eg states with both spin projections are unoccupied 
and tig-Bg hybridisation {t') lowers the energy of the FM 
as well as AFM configuration because in both cases only 
the bonding combination is occupied (see Fig. [TOl). The 
energy of the bonding state « — t'^/(£e^ — et2^) is, 
however, lower if the magnetisations of two Cr ions are 
parallel. In this case the partner with the same spin 
as the occupied tig state has lower energy Se^ = Ae_t vs. 



— Ae_f -I- A'li in the AFM case, so that 

2 2i'2A'^^ 
EpM - Eafm ^ 2JfmS ~ tt nrrr > C^) 

^e-H^e-t + ) 

and the effective coupling Jfm is ferromagnetic. It 
should be mentioned that these estimates for the effec- 
tive exchange couplings are in agreement with the well 
known Goodenough-Kanamori rules.— 

Thus, the strength and even the sign of Ji in ACriX^ 
are determined by the competition of two contributions 
of opposite signs, Jafm and Jfm, which depend on 
the strength of the tig-tig and tig-Cg hybridisation, re- 
spectively. In ZnCr204 Jafm dominates because of the 
strong direct hopping of dda type between the Cr tig 
states. This conclusion has been confirmed numerically 
by the results of test LSDA calculations in which the 
hybridisation between Cr d states was switched off. Ji 
derived from a fit to band energies calculated in this way 
is about 1 K, whereas for other J„ the values comparable 
to those given in Table |TT] were obtained. With the in- 
crease of the lattice constant in the row Zn^Cd^Hg the 
strength of the direct d-d hopping and, as a consequence, 
of Jafm rapidly decreases so that in IIgCr204 according 
to the LSDA results Jafm and Jfm almost exactly com- 
pensate each other. The dependence of Ji on the value 
of the parameter U in the LSDA-I-J7 calculations can be 
explained by weakening of Jafm caused by the increase 
of the denominator in Eq. ^ which in this case is given 
by A^i « 2 Jh + U. 

In y4Cr2S(e)4 spinels the distance between Cr nearest 
neighbors (see Table H]) is much larger than in ACr204 
and the contribution of the direct tig-tig hopping to Ji 
becomes less significant. For ZnCr2Se4, for instance, cal- 
culations with switchcd-off Cr d-d hybridisation give the 
value of FM Ji which is only 5% larger than the LSDA 
one. Thus, the indirect hybridisation via X p states be- 
tween the half-filled tig and empty Cg states split by the 
on-site exchange interaction provides the dominant con- 
tribution to FM Ji as was proposed by J. Goodenough 
in Ref. 18. A direct numerical proof is provided by test 
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calculations in which an external orbital-dependent po- 
tential Veg was added to Cr states; with its strength 



3Jh being equal to the LSDA exchange splitting 



of the Cg states. Then, the energies of Cr states with 
both spin projections are equal so that A'^-^ « and the 
nominator in the expression ([T]) for Jfm becomes zero. 
Ji — 15 K estimated from such calculations for ZnCr2Se4 
changes its sign to AFM and its absolute value is signif- 
icantly smaller than | Ji|=54 K obtained from the LSDA 
calculations. This indicates that the indirect t2g-X p-t2g 
hybridisation also contribute to Ji but it is much weaker 
than Jfm- Similar calculations for ZnCr204 result in 
substantial increase of the strength of AFM Ji which 
means the FM contribution is present also in oxides but 
J AFM wins due to direct d-d hybridisation. 

In order to explain the weird behavior of Ji in 
LSDA-|-f7 calculations one needs to analyze in more de- 
tails various parameters entering Eq. (O. Recalling that 
the orbital dependent part of the LSDA-I-J7 potential can 
be approximated as Vi — U'{l/2~ Ui) with U' = U — Jn 
one gets 



Ae-t 



3Jh + t^'(»T-T ^ ^l) ' 
AcF + U'{l-n^), 



(8) 
(9) 



where and ni are the occupation numbers of e^^ 
and states, respectively. Here we neglected the 
small (~0.05) deviation of the occupation numbers of the 
minority-spin t2g states from unity. Then, Eq. ([7]) can be 
written as 



2JfmS^ 



2t'2[3JH + C/'(nt -n^)] 
AcF + U'{l-n^) 
1 

icF + U'{l-ni)+3JH 



(10) 



If the occupations of the Cr were equal to zero Eq. 
([Toll would become 



2JfmS 



6t'2jH 



iAcF + U'){AcF + U' + 3Jii) 



(11) 



which would lead to fast decrease of Jfm with the in- 
crease of U'. However, since n| ~ 0.45 and ^ 0.25 
were calculated for ACr2S(e)4 the appearance of the term 
proportional to U'{ni — n^) in the nominator of Eq. pU]) 
slows down the decrease of Jfm- Also the terms propor- 
tional to n| and in the dominator of Eq. pT|) effec- 
tively decrease U', thus leading to the increase of Jfm 
as compared to purely ionic model with zero occupations 
of the Cg states. This is one of possible reasons why the 
FM contribution to Ji in overestimated in the LSDA-I-J7 
calculations. 

The calculated behavior of J2, J3, and J4 follows the 
predictions made by K. D wight and N. Menyuk in Refs. 
[19II20I on the base of the Goodenough-Kanamori rules of 
superexchange .^^1^'' Hybridisation paths Cr d-X p-X p- 
Cr d responsible for these interactions include an inter- 
mediate hopping between p states of two X ions. This re- 
duces the strength of the superexchange interactions but 



does not change the signs of t2g-t2g and t2g-eg contribu- 
tions to J„, which remain antiferro- and ferromagnetic, 
respectively. Indirect t2g-t2g and t2g-eg hoppings give 
comparable contributions to J2- Being of opposite signs, 
these contributions compensate each other. This explains 
the smallness of J2 and the variation of its sign in the 
ACT:2Xi spinels. Similar hybridisation paths between d 
states of 3-rd neighbors, with the dominant contribution 
of the indirect t2g-t2g hopping, lead to the appearance 
of AFM J3 and J3 . In the real structure the distance 
between the pair of X ions along the J3 path is smaller 
than along the J3 one; the difference being proportional 
to the deviation of the positional parameter x from 1/4. 
This may explain why the calculated values of J3 are 
somewhat larger than J3 . On the other hand, since the 
J3 pair of X ions belongs to the same AX^ tetrahedron, 
J3' is more sensitive to A s-X p hybridisation, which may 
be responsible for the increase of the J''^ / J'^ ratio in the 
Hg-containing spinels. Finally, indirect hybridisation be- 
tween t2g and Cg states gives dominant contribution to 
the superexchange between 4-th neighbors^^ which leads 
to FM J4. 



VII. SUMMARY AND CONCLUSIONS 

The LSDA(+J7) band structure calculations performed 
for some ACr2X4 spinels reproduce the change of the 
sign of the nearest-neighbor exchange interaction be- 
tween Cr spins from AFM in ^Cr204 compounds to FM 
in ACr2S(e)4 ones. It has been verified that the strength 
and the sign of Ji depend on relative strengths of two 
contributions of opposite signs, AFM and FM, which are 
due to Cr t2g-t2g and t2g-eg hybridisations, respectively. 
In ZnCr204 Jafm dominates because of the strong di- 
rect hopping between the Cr t2g states. As Zn is replaced 
by Cd or Hg, the strength of the direct hopping rapidly 
decreases with the increase of the Cr-Cr separation. As 
a result, Ji in the corresponding oxides is much weaker 
than in ZnCr204. In the ACr2S(e)4 spinels the FM con- 
tribution originating from the indirect hopping between 
half-filled t2g and empty Cg states, split by the on-site 
exchange interaction, becomes dominant. However, the 
magnitute of the FM contribution to Ji seems to be over- 
estimated by the calculations. The strongest among more 
distant exchange couplings are J3 and J3 to two sextets 
of 3-rd Cr neighbors. These interactions are always AFM 
and of comparable strengths. The couplings to 2-nd and 
4-th Cr neighbors are significantly weaker than Ji and J3. 
The obtained results are in accordance with the analy- 
sis of the exchange couplings in Cr spinels based on the 
Goodenough-Kanamori rules of superexchange and may 
be helpful for understanding the magnetic properties not 
only of Cr spinels but also Ti (MgTi204) and V (ZnV204, 
CdV204) ones, in which orbital degrees of freedom come 
into play^^iiSi^ 
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APPENDIX: 

Below explicite expressions for the difference Eb_{(1) ~ 
i?H(0), where -BhIq) is given by Eq. ^ and 

i;H(0) = 6Ji + 12J2 + 6J3 + 12J4, (A.l) 

are written as a functions of the angle (j) — qi: /2 for some 
directions of the wave vector q of spin spirals: 
q=(0,0,(?) {T~X): 

E{4>) = 4Ji(cos0- 1) +4J2(2cos(/) + cos2(/)- 3) 
+ 8J3(cos20- 1) 

+ 4J4(cos0 + cos30- 2); (A.2) 

q=(g,<Z,0) {T-K): 

E{(j)) = Ji(4cos0 + cos2(/i- 5) 

+ 2J2 (2 cos + cos 20 + 2 cos Si/) -5) 

+ 2J3(4cos20 + cos40- 5) 

+ 2J4(2 COS0 + cos 20 + 2 cos 30 + cos 40 — 6); 

(A.3) 



q_={q,q,q) (F-L): 



E{(j)) = 3Ji(cos20- 1) 

+ 3J2(2cos20 + cos40- 3) 

+ 6 Ja (cos 40- 1) 

+ 6J4(cos20 + cos40- 2); (A.4) 



q-(l,q,0): 



E{<j)) = 2Ji(cos0-3)-4J2(cos0 + 3) 

- 16J3 + 2J4(cos0 + cos30-6); (A.5) 



q=(i,i,g): 



E{(l)) = -6J1 - I2J2 

- 8J3(cos20+l)-12J4; (A.6) 



Here J3 = (J3 + J^)/2 is the average of the exchange 
coupling constants for two inequivalent groups of 3-rd Cr 
neighbors. 
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